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EFFECT  OF  THICKNESS  ON  THE  LATERAL  FORCE 
AND  YAWING  MOMENT  OF  A  SIDESLIPPING 
DELTA  WING  AT  SUPERSONIC  SPEEDS 
By  Kenneth  Margolis 

SUMMARY 


Calculations,  Based  on  linearized  supersonic -flow  theory  and  on  the 
additional  assumption  that  the  Mach  cones  yaw  with  the  wing,  are  made  to 
determine  the  effect  of  thickness  on  the  lateral  force  and  yawing  moment 
of  a  sideslipping  delta  wing  at  supersonic  speeds.  Generalized  equations 
are  derived  for  delta  wings  with  rhombic  profile  and  constant  thickness 
ratio.  Results  indicate  that,  for  wings  suitable  for  supersonic  flight, 
the  contribution  of  wing  thickness  to  the  stability  derivatives  Cyp 


and  Cnp  (rates  of  change  of  lateral-force  and  yawing-moment  coeffi¬ 
cients  with  angle  of  sideslip)  is  small  in  comparison  with  the  effects 
that  may  be  expected  from  a  vertical  tail. 


INTRODUCTION 


Theoretical  investigations  of  supersonic  stability  derivatives  of 
isolated  wings  have  been  confined  chiefly  to  the  treatment  of  mathemati¬ 
cally  thin  profiles  (flat  plates) .  Finite  thickness  will  have  certain 
effects,  some  of  which  can  be  estimated  by  consideration  of  small  pertur¬ 
bations.  The  purpose  of  the  present  paper  is  to  present  some  calculations 
showing  the  effect  of  thickness  on  the  lateral-force  and  yawing-moment 
derivatives  of  a  sideslipping  delta  wing  at  supersonic  speeds . 

On  the  basis  of  linearized  theory,  equations  are  developed  for  the 
stability  derivatives  Cyp  and  Cn  p  (derivatives  of  lateral  force  and 
yawing  moment  with  respect  to  sideslip)  as  functions  of  thickness  ratio, 
aspect  ratio,  and  Mach  number  for  delta  wings  at  zero  lift.  The  airfoil 
sections  parallel  to  the  plane  of  symmetry  are  rhombuses,  that  is, 
symmetrical  double  wedges  with  maximum  thickness  at  50  percent  chord,  and 
are  of  constant  thickness  ratio. 

In  the  analysis  the  Mach  cones  were  assumed  to  yaw  with  the  wing. 

This  additional  approximation  introduces  considerable  mathematical  simpli¬ 
fications  into  the  problem.  The  assumption  of  the  yawed  Mach  cones 
raises  a  question  as  to  the  quantitative  validity  of  the  results. 

Although  verification  of  these  quantitative  results  must  await  a  more 
exact  analysis,  the  qualitative  results  of  the  present  investigation 
are  believed  to  bring  out  the  salient  characteristics  and  signifi¬ 
cance  of  the  thickness  effects. 
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Generalized  curves  are  presented  to  permit  estimation  of  the 
derivatives  for  any  specific  delta  wing.  Series  of  curves  giving  varia¬ 
tions  of  the  derivatives  with  Mach  number  for  various  aspect  ratios  are 
also  presented. 


SYMBOLS 


General : 


V 

u 

w 

P 

Ap 

q 


Cartesian  coordinates 
flight  velocity 

incremental  disturbance  velocity  parallel  to  x-axis 
incremental  disturbance  velocity  parallel  to  z-axis 
density  of  air 

pressure  increment  caused  by  disturbance 


dynamic  pressure 


cp 

M 

B  =  \J  M2  -  1 
a 

c 


t 

A 

m 

b 

S 

A 


disturbance-velocity  potential 
Mach  number 


root  semichord,  measured  in  plane  of  wing  symmetry 

chord  length  at  spanwise  station  y,  measured  in  a 
plane  parallel  to  plane  of  wing  symmetry 

maximum  thickness  of  section  with  chord  c  at  spanwise 
station  y 

angle  of  sweep  of  leading  edge 
slope  of  wing  leading  edge  (cot  a) 
span  of  wing 
wing  area 


aspect  ratio 


cot  A 
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M-  =  cot  A^M2-1  =  ^  ^  M2  -  1 


P 

Y 

N 

Cy 

Cn 

R.P. 


angle  of  sideslip 
lateral  force 
yawing  moment 
lateral-force  coefficient 

yawing-moment  coefficient 

real  part 


|pV2S 


N 


Stability  derivatives : 


Subscripts: 

<  A 

yawing 

due 

<A 

yawing 

due 

arms 


moment  contribution  (taken  about  wing  apex) 
to  lateral  forces  multiplied  by  x  moment  arms 

moment  contribution  (taken  about  wing  apex) 
to  longitudinal  forces  multiplied  by  y  moment 


indicates  velocity  taken  in  flight  direction, 
except  when  used  for  stability-derivative 
notation  defined  previously 


ANALYSIS 


The  analysis  is  based  on  the  linearized  theory  of  supersonic  flow 
and  the  derivation  of  equations  and  application  of  results  are  therefore 
subject  to  the  usual  limitations  and  restrictions  of  the  theory.  In  order 
to  minimize  the  mathematical  complexities,  all  pertinent  equations  are 
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derived  specifically  for  a  delta  plan  form  with  rhombic  profile.  The 
method  is  applicable,  however,  to  wings  of  arbitrary  plan  form  and  profile 

Consider  the  delta  wing  shown  in  figure  l(a) .  The  wing  is  moving 
in  a  direction  parallel  to  the  axis  of  symmetry,  that  is,  at  zero  lift 
and  zero  sideslip.  The  incremental  disturbance  velocity  parallel  to 
the  x-axis,  u,  is  readily  obtained  by  methods  outlined  in  references  1 
to  3.  The  equation  is 


(1) 


where  —  and  y  for  the  ^omrarA.  wedge  and  rearward  wedge,  respectively 
are  found  as  follows: 


It  should  be  noted  that  equation  (l)  is  applicable  to  other  Mach  line 
configurations  as  well  as  for  that  shown  in  figure  1. 

Wow  consider  the  same  wing  sideslipping  to  the  right  at  a  small 
angle  p.  (See  fig.  1(h).)  The  Mach  cone  is  assumed  to  retain  the  same 
position  relative  to  the  wing  for  both  the  straight-flight  and  the  side¬ 
slipping  condition.  This  assumption  follows  from  consideration  of  the 
discussion  given  in  reference  k.  The  expression  for  u  for  the  side¬ 
slipping  wing  is  still  given  by  equation  (l)  except  that  other  values  are 

required  for  —  and  y.  For  small  angles  of  sideslip  these  values  are 
approximately  as  follows: 

For  the  right  half -wing. 
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and  for  the  left  half -wing. 


T  =  ~(l  -  P  tan  A) 


Y  3  ~  P  tan  A) 


For  the  sideslipping  condition,  equation  (l)  "becomes 


8$  =  b.P.  2S£  i 
ox  «  c 


2  2 

1  -  m  B 


cosh"1  £  liD 
B|y  -  mx| 


;.l  ~  M°°A)  oosh-l  X  +  mBgT 

-ll  -m2B2  +  H 


4m2B2 


cosh-1  a  -  2mE2; 
B|y  -  2m(x-a 


ljm2B2 


-1  x  -  a  +  2mB  ,y 
B|  y  +  2m(x-a) 


The  disturbance  velocity  in  the  flight  direction  Up  (positive  forward) 
may  be  expressed  in  terms  of  components  along  and  perpendicular  to  the  axis 

of  wing  symmetry,  ^2  and  42 .  The  relation  is 
dx  dy 


The  increment  in  pressure  caused  by  the  disturbance  is  obtained 

from 

Ap  =  pVup  (6) 

which  expresses  Bernoulli's  law  with  the  approximation  of  small  disturbances. 

From  equations  (5)  and  (6),  the  pressure  coefficient  ^  for  one 
surface  of  the  sideslipping  delta  wing  is  found  to  be 


Ap  =  2/42  -  r  42 

q  V\Bx  By. 
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The  expression  for  is  given  in  equation  (4) }  the  term  may 

ox  oy 

“be  obtained  in  the  following  manner: 

82  =  A  82  dx  (8 

by  by  8x 


The  expression  for  is  therefore 


=  b.P.  i  V  i  1  -  |  tanA  Cosh_1  ■g.+  .Pgl 


*  0  2  \  1  -  m2B2 


B  IE  + 


t  P  tan.  A  cosh- i  x.  +  p  tan  A  (cosh 


2  Ml  -  m2B^ 


B  mx  -  y| 


-1  x 


cosh 


2  -  (3  tan.  A  cosh-l  x  -  a  +  2aiB 


2  \  1  -  4nrB‘ 


B|2m(x  -  a)  +  yl 


2  +  3  tan  A  _ -1  *  -  a  -  2nTJ  (> 

r  ■  ~  77  B|a.(x  - «)  -  d 

2  Ml  -  tort2 

Performing  the  operations  indicated  in  equation  (j)  and  neglecting 
all  terms  of  order  p2  and  higher  gives  the  following  simplified 
expression  for  the  pressure  coefficient  of  the  sideslipping  delta  wing: 


Ag  _  E  p  2  t  m  +  2p 

q  ”  '  ' 11  c  \TT^ 


cosh-1  ^y,  +  m  Z  2p  cosh-1  x,+  Bg5 


Bly  - 


2  2 

1  -  m B 


By  +  mx 


i(ia  +  cosh-1  x a  ;4b^ 


rw  Bly  ‘  2m(X_a)l  \|l  -  4m2 B2 


-  11  cosh"1  x,-  a  + 


B|y  +  2m(x-a)| 
(10-)' 


The  lateral- force  and  yawing-moment  coefficients  may  then  he  obtained 
by  appropriate  integrations  over  the  entire  wing  area,  thus , 


Cv  = 


lpT2S 


dx  dy 

q  cLy 
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and 


:  =  l 

n  ipV2Sb 
2 


2 

Sb 


q  <V  J  J  q  d2 


y  dx  dy 


(12) 


where,  for  convenience,  the  moments  are  taken  about  the  wing  apex  and  are 
considered  positive  for  clockwise  rotation.  The  factor  of  2  accounts  for 
the  equal  additive  contribution  of  the  upper  and  lower  surfaces.  Informa¬ 
tion  pertinent  to  the  limits  of  integration  is  shown  in  figure  2.  The 
integral  expressions  for  CY  are  presented  in  appendix  A  and  those  for  Cn 

in  appendixes  B  and  C .  Appendix  B  presents  the  first  term  on  the  right- 

hand  side  of  equation  (12)  which  is  designated  (Cn)  j  appendix  C  presents 

the  second  term,  which  is  designated  (Cn)  .  x 

a 


The  indicated,  integrations  were  performed  and  the  resulting  expressions 
for  the  stability  derivative  Cyp  and  for  the  derivatives  {^n^)x 

and  are  given  as  functions  of  p.  in  table  I* 


For  moments  taken  about  an  arbitrary  point  Xq 
may  readily  be  derived  as 


an  expression  for 


%  -  (s)x  +  (%)y  +  I  1?  %  (!3) 

Values  of  the  stability  derivative  Cn^  are  then  obtained  by  adding  the 
previously  derived  right-hand  members  of  equation  (13) . 


RESULTS  AND  DISCUSSION 


It  was  pointed  out  in  the  introduction  that  the  assumption  of  yawed 
Mach  cones  raises  a  question  as  to  the  quantitative  validity  of  the 
results.  The . small  triangular  regions  of  the  wing  between  the  assumed 
(yawe  )  position  of  the  Mach  cones  and  the  correct  (unyawed)  positior 

have  an  area  proportional  to  the  sideslip  angle  0.  It  appears  that  the 

basic  pressure  for  0=0  acting  on  these  small  triangular  regions  may 
make  a  contribution  to  the  stability  derivatives  and  Cnp  that  is 

of  the  same  order  of  magnitude  as  that  produced  by  the  increment  in 
pressure  due  to  0  acting  over  the  entire  wing.  The  quantitative  results 
ob  tamed  ^herein  must  therefore  be  considered  as  rough  approximations. 

The  qualitative  results,  however,  which  represent  the  more  important 
consideration  of  the  present  investigation,  are  believed  to  bring  out  the 
salient  characteristics  and  significance  of  the  thickness  effects. 
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Stability  derivative  Cy  A  generalized  curve  showing  the  variation 


of  a  — with  |i  where  n  =  -  \|m  -  1  is  presented  in  figure  3  •  It 

(t/c)2  4 

should  he  noted  that  the  two  discontinuities  in  slope  evident  in  figure  3 
(also  in  figs.  4  to  6)  occur  when  the  Mach  lines  are  parallel  to  the 
maximum- thickness  line  and  to  the  leading  edge.  For  any  given  values  of 
aspect  ratio,  Mach  number,  and  thickness  ratio,  the  stability  deriva¬ 
tive  CY(3  is  readily  obtained  by  use  of  this  single  curve.  Variations 

Cy 

of  P  -  with  Mach  number  for  various  aspect  ratios  are  presented  in 
(t/c)2 

figure  k*  The  contribution  of  wing  thickness  to  the  value  of  Cy^  for 

suitable  supersonic  wings  (10  percent  thick  or  less)  is  found  to  Ido  small 
in  comparison  with  the  value  of  Cy^  that  may  be  expected  from  a  vertical 


The  contribution  of  wing  thickness  to  the  value  of  Cy^  for 


tail. 


Stability  derivative 


(t/c) 


(t/c) 


-  Generalized  curves  showing  the  variations 


with  [i  are  shown  in  figures  5(a)  and  5(b), 


respectively.  For  any  given  set  of  parameters,  the  stability  deriva¬ 
tive  Cnp  is  obtained  by  proper  use  of  figures  k  and  5  and  equation  (l3 )  • 


Variations  of  - ft-  (moments  taken  about  the  --chord  point)  with  Mach 

(t/c)2  3 

number  for  various  aspect  ratios  are  presented  in  figure  6.  As  with,  Cy^ 
the  contribution  of  wing  thickness  to  Cnp  for  wings  of  small  thickness 

is  seen  to  be  small  compared  with  the  expected  contribution  of  a  vertical 
tail,  although  for  some  airplane  designs  the  value  of  Cnp  due  to  thick¬ 
ness  may  be  appreciable  when  compared  with  the  value  of  Cn 


complete  airplane . 


for  the 
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CONCLUDING-  REMARKS 


On  the  basis  of  linearized  theory,  calculations  have  been  made  to 
determine  the  effect  of  thickness  on  the  lateral  force  and  yawing  moment  of 
a  sideslipping  delta  wing  at  supersonic  speeds.  Results  indicate  that,  for 
wings  suitable  for  supersonic  flight,  the  contribution  of  wing  thickness  to 
the  values  of  Cyp  and  Cnp  is  small  in  comparison  with  the  effects  that 


may  be  expected  from  a  vertical  tail,  although  for  some  airplane  designs 
the  value  of  Cf  due  to  wing  thickness  may  be  appreciable  when  compared 
np 


to  the  value  of  C~ 

np 


for  the  complete  airplane . 


Langley  Aeronautical  Laboratory 

National  Advisory  Committee  for  Aeronautics 

Langley  Air  Force  Base,  Va.,  October  6,  1948 
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Case  III 

4p(t/c)2 


(n  >  1) ; 


ni  _ 


fy/m 


rt  cbc  dy 


'2am 


y+2am 

2m" 


2am  „  y+2am 


2mB-l 


2m 


2am  J  y/ffl 
2mB-l 


it  dx  dy  + 


yB 


D  (k  dy 


D  = 


cos 


-1 


p 

x  -  mB  y 
B|y  -  mx| 


cos 


-1 


2 

x  +  mB  y 
B|y  +  mx| 


where 
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APPENDIX  B 

INTEGRAL  EXPRESSIONS  FOR  (Cn)x 


Integral  expressions  for  the  yawing- moment  contribution  (Cnj. 
for  the  three  cases  shown  in  figure  2  are  as  follows: 


Case  I  (0  <  [x  <  0: 

jta^m^ 


y+2am 
2  am  /  2m 


p(t/c) 


^  r2am  r  rc 

Ml  -  m2B2  h  h/m 


Ax  dx  dy 


am 


y+2am 


y+2am 


1-mB  /  2m 


Cx  dx  dy  + 


By+a 


f 2am  r  2m 

am  J  y/m 

1-mB 


Cx  dx 


where 


and 


2  2 
n  -1  x  -  mB  y  .1  x  +  mB  y 

A  =  cosh  -1  — - ■Lr  -  cosh  x  — — 7 — 

B  y  -  mx  B  y  +  mx 


Q  O 

-1  x  -  a  -  2mB  y  ,  -1  x  -  a  +  2mB  y 

C  =  C°Sh  B|y  -  2m(x  -  a)  |  °°S  Bjy  +  2m(x  -  a)  | 


Case  II  (J  <  ix  ^  l): 


00  i  A2am  f  ^ 

JtabmP  fn  \  -  -  _ ± _  /  /  Ax  dx  dy 

J  0  V  y 


MT  m2B2 


p(t/c)2 

where'  A  is  the  function  defined.  for  case  I 


y+2am 

r 
c 

y/m 
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Case  III  (n  >1) 
P(t/c): 


(°n)v  '  - 


\|m2B2  - 


2am 

2mB-l 


yB 


y/m 


rtx,dx  dy 


p  2am 

2m 

2mB-l 

2m  \ 

/  1 

jtx  dx  dy  + 

Dx  dx  dy] 

1  2am  J 

y/m  * 

0  J 

yB  / 

2mB-l 


where 


D  = 


cos 


x  -  mBgy 
B|y  -  mx 


cos 


-1 


x  +  mBgy 
B|y  +  mxl 


11*. 
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INTEGRAL  EXPRESSIONS  FOR  (Cn)i 


Integral  expressions  for  the  yawing-moment  contribution 
for  the  cases  shown  in  figure  2  are  as  follows: 


Case  I  (0  u  < 


-*.ais.  -  /c  \  = 

P(t/c)2  ~ 


1  -  m2B2 


y+2am 

'  2am  C  2m 
0  J  y  /m 


2am  P2a 


Ay  dx  dy  - 


Ay  dx  dy 


y+2am 

“ST” 


1  -  4m2®2 


(r  2am  ^  2a 

J  0  J  y+2am 
2m 


Cy  dx  dy 


y+2am 

^  2am  f  2m 

am  J  y/m 
1-mB 


Cy  dx  dy  - 


am  y+2am 
1-mB  f  2m 

1  J  By+a 


Cy  dx  dy 


where 


A  .  oosh-1 

B |y  -  mx 


-1  x  +  mB  y 

L  - 

B  y  +  mx 


C  =  cosh 


-1  x  -  a  -  2mB^y  x  -  a  +  2mB^y 

"*■  _ _ _ _ _ _ _ _ _ _  _  /1  X  -  -  -  - 


B|iy  -  2m(x  -  a)j 


B|y  +  2m(x  -  a)| 
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where  A  is  the  function  defined  for  case  I  and 


E 


cos 


-1 


x  -  a  -  2mB2y 
By-  2m (x  -  a)| 


cos 


-1 


x  -  a  +  2mB2y 
B| y  +  2m(x  -  a)| 
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FOR  LATERAL- FORCE  AND  YAWING-MOMENT  DERIVATIVES 


y  -  value 


Figure  3.-  Generalized  curve  presenting  the  variation  of  lateral -force 
derivative  with  a  parameter  relating  aspect  ratio  and  Mach  number* 


Variation  of  lateral -force  derivative  with  Mach  number  for  various  aspect  ratios 


